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ABSTRACT 

We present three-dimensional hybrid simulations of coUisionless shocks that propagate parallel to 
the background magnetic field to study the acceleration of protons that forms a high-energy tail on 
the distribution. We focus on the initial acceleration of thermal protons and compare it with results 
from one-dimensional simulations. We find that for both one- and three-dimensional simulations, 
particles that end up in the high-energy tail of the distribution later in the simulation gained their 
initial energy right at the shock. This confirms previous results but is the first to demonstrate this 
using fully three-dimensional fields. The result is not consistent with the "thermal leakage" model. 
We also show that the gyrocenters of protons in the three-dimensional simulation can move away from 
the magnetic field lines on which they started, which is not the case in one-dimensional simulations. 
Our study clarifies the injection problem for diffusive shock acceleration. 

Subject headings: acceleration of particles - cosmic rays - shock waves - turbulence 



1. INTRODUCTION 

CoUisionless shocks and the associated accelera- 
tion of charged particles are among the most im- 
portant processes in spac e physics and astrophysics 
i|Blandford fc Eichleii HQSTI). The theory of diffusive 
shock acceleration (DSA: iKrvmskvl Il977t lAxford et al.l 
[l977t [Belli 119781: IBlandford fc Ostriked llQTST describes 
quantitatively the acceleration process in the vicinity of 
shock waves. It predicts a power-law energy spectrum 
downstream of the shock, which is a common character- 
istic of the energetic charged particles observed in space. 
However, DSA does not address the physical processes 
involved with the acceleration of thermal and/or low en- 
ergy particles which, presumably, constitute the source 
of the high-energy particles. This is known as the "injec- 
tion problem" which has received a lot of recent attention 
but, as yet, t here is not a commo n consensus on its res- 
olution fe.g.. lKirk fc Dendvl[200TI) . 

CoUisionless shocks propagating in magnetized plas- 
mas are usually divided into two classes based on the an- 
gle, 9Bm between the incident magnetic field vector and 
shock normal vector. Quasi-parallel shocks are those for 
which 0° < Obtl < 45° and quasi-perpendicular shocks 
are for 45° < Bbu < 90°. For quasi-parallel shocks, 
there have been several mechanisms proposed in order 
to solve the injection problem. It is usually thought that 
for parallel shocks Alfven waves excited by the stream- 
ing of protons can scatter ions in pitch angle. The par- 
ticles can be acce lerated by DS A when their anisotropy 
is small enough. lEUisonl (|1981[ ) first proposed a Monte 
Carlo model for DSA that includes the injection process, 
where the injected particles originate from the shock- 
heated ions downstream of the shock that then leak back 
upstream of the shock. Subsequently, the particles are 
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"injected" into the standard DSA process of scattering 
back and forth across the shock. The model assumes that 
all the particles are scattered by magnetic fluctuations 
and ignores the details of the shock structure. It has been 
used t o fit and compare with observatio ns at Earth's bow 
shock (|Ellisonlll985l : lEUison et al.lll990( ). Similar models 
have also been discusse d by other authors (e.g.. iMalkovl 
Il998t iKang et al.ll2012[ ). This is usually referred to as 
the "thermal leakage" model. 

However, a number of authors have found a dif- 
ferent scenario for the initial energization at par- 
allel s hocks based on self-consistent hybrid simula- 
tions (lOuestl [l98l iScholer k Terasawal fl990t iScholeii 



19901 IKucharek fc Scholcr" 1 19911: IGiacalone et al.l lT992t 
Sugivama fc Ter asawa 1999: iSu et al.ll2012D . Thev find 



that the accelerated ions originate from the shock tran- 
sition layer rather than the downstream region. Some 
initially thermal ions are accelerated to high energies 
as they execute cycloidal motion within the electric and 
magnetic fields in the shock layer. Although the average 
incident magnetic field is parallel to the shock normal, 
as the enhanced upstream magnetic fluctuations steepen 
and convect through the shock layer, the angle between 
the incident magnetic field and the shock normal right 
at the shock front can significantly deviate from the av- 
erage value and ev en be close to a quasi- perpendicular 
shock locally (e.g., iWilson HI et al.ll2013() . A particle 
can gain the first amount of energy b y moving against the 
electric field within the shoc k layer ( Kucharek fc Scholerl 
119911: IGiacalone et al.lll992f). It has been clearly shown 
bv lKucharek fc Scholed ([1991 that most of the acceler- 
ated particles are reflected by the shock and gain the first 
increment o f energy at the sho ck layer itself. It is note- 
worthy that iLvu fc KanI ()199GD also used hybrid simula- 
tions and found that the leakage of protons dominated 
the accelerated particles. We suggest that the reason 
they obtained a different result from other authors may 
be due to the method they used to drive shocks. In their 
simulations, the upstream and downstream plasmas are 
initially connected smoothly by assuming a hyperbolic 
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tangent function for the shock parameters with a thick- 
ness of several ion inertial lengths. The magnetic fluctu- 
ations that are important to reflect ions at shock front 
arc ignored at the beginning of the simulation. Since the 
simulation only lasts for 85r2~/, where fid is the gyrofre- 
quency of protons in the asymptotic upstream magnetic 
field, it is difficult to determine if the leakage of protons 
can dominate the accelerated population of particles af- 
ter the shock is fully developed. There is observational 
evidence in support of the idea that particles originate in 
the shock layer. For example, signatures of ion reflection 
at quasi- parallel portion of Earth's bow shock has been 
reported (jGosling et al.lll982L [19891) . It is worth noting 
that although the initial acceleration in these two pro- 
cesses are quite different, these two models are closely 
connected and it is difficult to distinguish between them 
in observations. The Monte Carlo model has been ex- 
plicitly compared with the results of the hybrid simula- 
tions. The comparison genera hy has a good agre ement 
but some obvious differ ences (lEllison et al. I1993D . The 
readers are referred to ([Burgess et al.l 120121 ) for a com- 
plete review of these two models. 

The injection process at quasi-perpendicular shocks 
is thought to be more complicated. The required pre- 
acceleration may be achi eved by processes such a s shock 
drift acceleration (e.g., [Armstrong et al.l 119851: iDeckerl 
1988") and shock surfing acceleration ([Lee et al .1 119961 : 
Zank ct al. 199fi). In shock drift acceleration, charged 
particles drift because of the gradient in the magnetic 
field at the shock front. The direction of the drift is 
in the same direction as the motional electric field vec- 
tor E = —V X B/c, and the particles gain energy dur- 
ing this drift motion. In shock surfing acceleration, it is 
thought that the cross shock potential electric field is a 
barrier to some particles incident on the shock leading to 
their reflection at the shock and a drift against the mo- 
tional electric field leading to the particle acceleration. 
Recent progress has been made to distinguish the rela- 
tive importance for these two processes. For instance, 
it has been found that in order for shock surfing accel- 
eration to be efficient, the thickness of the shock layer 
has to be very thin - on the s cale of an electron inertial 
length ([Lipatov fc Zanklll999D . which is not consistent 
with the shocks observed i n space and in nu merical sim- 
ulations ([Bale et al.ll2003[: iLerov et al.lll982[ ). Moreover, 
it has been shown that the shock thickness has to be 
fairly large compared to electron gyroradii to be consis- 
tent with the observa tion of electron heating at shocks 
([Lembege et al.|[2004[ ). More detailed studies using full 
particle sim ulations and hybri d si mulations have b een 
presented by 'Yang ct al." and [Wu et ahl (|2009t) . 

Giacalone & Jokipii (199^ have demonstrated that 
large-scale magnetic turbulence which is present in the 
background plasma can efficiently lower the injection 
threshold at perpendicular shocks by increasing the 
transport of charged particles normal to magnetic field. 
Recent numerical simulations for the acceleration of 
charged particles (both ions and electrons) in the ex- 
istence of large-scale magnetic fluctuations show very 
efficient acceleration, which indicates that t here is no 
injecti on problem at perpendicular sh ocks ([Giacalond 
[2005al[bl : IGuo fc GiacaM[20Tol [20T2a[ [bl). It is worth 
noting that observationally there seems no injection 



problem at quasi-perpendicular shocks. For instance, 
the Voyager spacecraft clearly show energetic particles 
(at least several M eV) are accelerated at th e solar wind 
termination shock ([Decker et al.l [20051 12008| ). Strong in- 
terplanetary shocks observed by ACE and Wind at 1 AU 
(mostly quasi-perpendicular shocks ] can accelerate pro- 
tons to energies more than 50 keV ([Giacaloiii[2012|) . 

It is important to point out that, in order to ac- 
curately model the acceleration of charged particles at 
shocks, one has to consider fully three-dimensional elec- 
tromagnetic flelds. This is because in a field that has 
at least one ignorable coordinate, the motions of parti- 
cles are artificially restricted in the sense that any in- 
dividual charged particles must remain within one gy- 
roradius o f the magnetic field line on which it began 
its motion (Jokipii et al.lll993[ : [Giacalone fc Jokipiil[l994[ : 
[Jones et al.l[1998| ). Many prev ous hybrid simulations 
were performed using only one or two spatial dimensions 
(Scholer et al., 1993), and, thus, are subject to this con- 
straint. The effect of cross-field diffusion must be in- 
cluded and requires fully three-dimensional simulations 
([Giacalone fc Ellison[[2000[ ). 

In this study, we employ a three-dimensional self- 
consistent hybrid simulation to study the initial particle 
energization at a parallel shock. We find that all parti- 
cles accelerated to high energies gain the first amount of 
energy right at the shock front. In the three-dimensional 
simulation, the accelerated protons can move off their 
original field lines. Nevertheless, the initial acceleration 
process is consistent with that is found in previous one- 
dimensional and two-dimensional hybrid simulations. 

2. NUMERICAL METHOD 

We perform three-dimensional hybrid simulations of 
parallel shocks and compare the results with that of one- 
dimensional hybrid simulations. Although our simula- 
tions resolve the full proton distribution function, we 
are particularly interested in those particles that are ac- 
celer ated to high en ergies. In the hybrid simulations 
(e.g.. [Winske fc Que st 1988), the ions are treated kinet- 
ically and thermal electrons are treated as a massless 
fluid. This approach is well suited to resolve ion-scale 
plasma physics, which is crucial to describe supercrit- 
ical coUisionless shocks and particle acceleration from 
thermal to suprathermal energies. We have improved 
the efficiency of our one-dimensional, two-dimensional, 
and thre e-dimensio nal hybrid simulation models (Section 
3.3.1 in [GMl[20T2[ ). The new codes have been imple- 
mented and tested on the NASA's Pleiades supercom- 
puter using a few thousand CPU cores. For the three- 
dimensional simulations, we consider a three-dimensional 
Cartesian grid (x, y, z). All the physical vectors such as 
the positions and velocities of protons, and electric and 
magnetic flelds E and B have components in three direc- 
tions and also spatially depend on x,y, and z. For the 
one-dimensional simulations, the physical vectors have 
components in three directions but only depend on x. A 
shock is produced by using the so-called piston method, 
in which the plasma is injected continuously from one 
end of the simulation box {x = 0, in our case), and re- 
flected elastically at the other end (x = L^)- The right 
boundary is assumed to be a perfectly conducting bar- 
rier. The pileup of density and magnetic fleld creates 
a shock propagating in the —x direction. In the three- 
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dimensional simulation, the boundary conditions of the 
electromagnetic fields in the y and z direction are peri- 
odic, and the particles that move out of one end of the 
simulation domain in y or z direction will re-enter the 
domain from the other end. 

We examine two different simulation cases (one- 
dimensional and three-dimensional) with similar param- 
eters. In both cases, the Mach number of the flow in the 
simulation frame is M^o = Vx/Vao =4.0, where Vao is 
the upstream Alfven speed. The electron and ion plasma 
betas are /3e = 1.0 and (ii = 0.5, respectively. The grid 
size is Aa; x Ay x Az = O.Sc/wpi x 0.5c/wpi x O.Sc/wpi 
for the three-dimensional case and Ax = 0.5c/wpi for the 
one-dimensional simulation case, where c/wpj is the ion 
incrtial length, and c and ojpi are the light speed and 
proton plasma frequency, respectively. The time step is 
taken to be Vld^t — 0.01, where Vld is the proton gyro- 
frequency. The ratio between light speed and upstream 
Alfven speed is c/vao = 6000.0, and the anomalous resis- 
tivity is 77 = 1 X 10~^47ra;~"'^. The initial spatially uniform 
thermal ion distribution is represented using 25 particles 
per cell in the three-dimensional case and 200 particles 
per cell in the one-dimensional case. We have repeated 
the simulations with different number of particles per cell 
and found it docs not change the results. Initially the av- 
erage magnetic field is assumed to be Bq ~ Bqx^ i.e., the 
average shock-normal angle is zero. The spatial size in 
the X direction is taken to be 300c/cjpi for both of the 
cases. For the three-dimensional case, the lengths of the 
simulation box in the y and z directions are taken to be 
Ly = Lz = 'jQc/u Jm- Pre-existing magn etic fluctuations 
are not included (jGiacalone et al.lll992| ). 



3. SIMULATION RESULTS 

Figure [T] shows results from the three-dimensional sim- 
ulation. Shown are representations of various quantities 
in the x-y plane at z = 20 c/tUpi at Q,cit ~ 120.0. The 
figure shows the color-coded contours of (a) the magni- 
tude of the magnetic field B/Bq, (b) the x-componcnt of 
the bulk proton velocity Vx/Vao, (c) the plasma density 
n/no and (d) the density of accelerated particles he /no 
with energies of 3Ei < E < 5Ei, where Ei represents 
the upstream proton ram energy mpUl/2 in the shock 
frame, and Ui is taken to be 5.3Vao, which is the up- 
stream plasma flow speed measured in the frame of the 
shock. From these plots we can see that the shock front 
is approximately located between x = 140 - 150c/ LUpi 
where the ion density and flow speed have clear, abrupt 
transitions. The magnetic field and fluid velocity in the 
shock region show strong fluctuations. The fluctuations 
also modify the spatial distribution of plasma density 
and the density of accelerated particles. The accelerated 
protons are concentrated around the shock layer where 
the magnetic field increases, indicating a strong localized 
energization in this region. 

In Figure [2] and Figure [3] we compare the x-profiles of 
the three-dimensional hybrid simulation with that of the 
one-dimensional hybrid simulation. The black solid lines 
in Figure [2] show (a) the plasma density n/riQ, (b) the 
magnitude of magnetic held B/Bq, and (c) flow veloc- 
ity in the x direction Vx/Vao along y ~ 20 c/ujd and 
z = 20 c/u)ci at Vlcit = 120.0. The red dashed lines show 
the x-profiles from the one-dimensional simulation at 
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Fig. 1. — Two-dimensional representation of the three- 
dimensional hybrid simulation in the 2 = 20 c/u)pi plane at 
Q,cit = 120.0: (a) the magnitude of the magnetic field B/Bq, (b) 
the x-component of the bulk proton velocity Vx/Vao^ (c) the den- 
sity of the plasma flow n/riQ, and (d) the density of the accelerated 
particles with energy SSi < E < 5Ei , where Ei represents the up- 
stream proton ram energy mpU^/2 in the shock frame, and Ui is 
taken to be S.SV^o- 



Vlcit — 110.0. The black solid lines in Figure [3] show the 
x-profiles from the three-dimensional simulation along 
the same line with that in Figure [2J but for the magnetic 
field components in Bx/Bq, By/ Bq, and B^/Bq, respec- 
tively. The red dashed lines represent the results from 
the one-dimensional simulation at flcit = 110.0. For the 
onc-dimensional simulation, the location of the shock in 
the x direction is about x ~ 132c/ajci and the upstream 
inflow speed in the shock frame is about Ui = 5.5Vao, 
slightly larger than that in the three-dimensional sim- 
ulation. For the three-dimensional case, the magnetic 
field fiuctuations are mostly transverse to the initial mag- 
netic field in the x direction. The fluctuating components 
5By and SB^ contain about 90% of the total energy of 
the magnetic fluctuations. The fluctuations are circu- 
larly right-handed polarized, and are convected toward 
the shock, indicating that they are excited by reflected 
ions that flow upstream. Close to the shock, the fluctu- 
ations are compressed and amplifled. The downstream 
magnetic field consists of large-amplitude fluctuations. 
The results are consistent with previous analytical the- 
ories aiid_one-dimcnsional numerical simulations (e.g., 
iQuestI [l988f ). The results from three-dimensional sim- 
ulations are qualitatively consistent with those from the 
onc-dimensional simulation. 
Now we discuss the acceleration of protons at the shock 
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Fig. 2. — One-dimensional x-profilcs from the three-dimensional 
simulation at j/ = 20c/u)pi and z = 2f)c/ujpi (black solid lines) 
and the one-dimensional simulation (red dashed lines): (a) plasma 
number density n/no, (b) the magnitude of the magnetic field 
B/Bq, and (c) the x-component of the bulk proton velocity 
V-^/Vao- 



fron t. iJokipii et all (I1993D . iGiacalone fc Jokipiil ()1994[ ) 
and lJones et al.l (|1998^ ~ have demonstrated that in simu- 
lation models with at least one ignorable coordinate, the 
guiding center of an individual charged particle is con- 
fined to within one gyroradius of the magnetic field line 
on which it began its motion. In order to overcome this 
artificial constraint and to accurately model the acceler- 
ation process, it is important to consider the motions of 
charged particles in three-dimensional electric and mag- 
netic fields. Since we perform three-dimensional hybrid 
simulations in this study, the restriction on the motions 
of charged particles is removed. 

Figure |4] shows the trajectory of a representative par- 
ticle in the one-dimensional simulation. The physical 
quantities for each panel are: (a) the energy versus the 
location in the x direction away from the shock front in 
the shock frame Xsh (the location of the shock front is 
assumed to be Xsh = 0), (b) the energy versus the time 
in the simulation, (c) the simulation time versus Xsh , and 
(d) the location in the z direction z ~ J v^dt -\- zq versus 
the location in the y direction y = J v.ydt + yo, where j/o 
and Zq are assumed to be zero, respectively. The marker 
'1' notes the starting point and '2' notes the end point 
of the trajectory. One can clearly see that the initial ac- 
celeration occurs right at the shock front Xsh = 0. The 
acceleration mechanism is due to the reflection at the 
shock, and the particle gains energy in the electric field 
E = —V X B/c since locally the angle between magnetic 
field vector and shock normal vector can be quite large. 
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Fig. 3. — Similar to Figure [J] but for x-profiles of (a) the x- 
component of the magnetic field B^/Bq, (b) the y-componcnt, 
By/Bo, and (c) the z component, Bz/Bq. 
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Fig. 4. — The trajectory of a representative particle accelerated at 
the shock layer from the one-dimensional simulation. The physical 
quantities for each panel are: (a) the energy versus the distance 
away from the shock layer, (b) the energy versus the time in the 
simulation, (c) the simulation time versus the distance away from 
the shock, and (d) the location along the z direction z = J v^dt+zo 
versus the location along the y direction y = J Vydt + yo, assuming 
the initial location is at j/o = and zq = 0. 



As shown in this figure, the particle "rides" along the 
shock front and is accelerated for about 15 gyroperiods. 
The resulting energy gain is about 10 times of the plasma 
ram ener gy. Th is process h as been discusse d by previous 
authors (lQuestlil98& .Scholcr fc Terasawal,1990; .Schole^ 
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Fig. 5. — The trajectory of a representative particle accelerated 
at the shock layer from the three-dimensional simulation. The 
plotted physical quantities are similar to Figure |4] 
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Fig. 6. — The locations of 50 accelerated particles in the x- 
direction as a function of simulation time. The particles are tracked 
until they gain energies several times the plasma ram energy. It 
shows that all the accelerated particles gain initial energies at the 
shock layer. 



[19901 IKucharek fc Scholed [T99lt IGiacalone etHI [1991) 
and our results are consistent with the previous work. 
Since we use the one-dimensional hybrid simulation, the 
motion of the particle is restricted on its original field 
line and its guiding center to within one gyroradius of 
the initial field line. 

In Figure \5\ we show the trajectory of an accelerated 
particle in the three-dimensional simulation. The plotted 
physical quantities are similar to Figure|4l Since the elec- 
tric and magnetic fields are fully three-dimensional, this 
particle (and all others) is not restricted to its original 
field line as in the one-dimensional simulation. This can 
be seen from Figure [5] (d), in which the guiding center 
of the charged particle drifts in the y-z plane away from 
its original field line during the acceleration. We have 
examined the trajectory of this particle until the end of 
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Fig. 7. — The displacements of particles from their original po- 
sitions in the z-direction versus the location of particles in the 
x-direction. The upper panel shows the results from the one- 
dimensional simulation. The bottom panel shows the results 
from the test-particles that are initially on the same field line 
y = 20c/ajpi and z = 20c/LUpi in the three-dimensional simulation. 



the simulation and found that its guiding center does not 
move back to within one gyroradius of the position of the 
field line that the particle started its gyromotion. The 
original field line is not likely far away (in the y-z plane) 
from where it was originally. In order for this to happen, 
there would have to be significant flows in the y and z 
direction upstream of the shock. 

However, even in the three-dimensional simulation we 
find that the acceleration process is quite similar to what 
is found in the one-dimensional simulation. The particles 
still gain their first amount of energy at the shock front. 
We have examined several hundred accelerated particles 
in the three-dimensional simulation and found that they 
all get the initial acceleration right at the shock layer. 
In Figure [Sj we randomly selected 50 accelerated parti- 
cles in the three-dimensional simulation and show their 
positions in the x direction as a function of time. The 
particles arc tracked until they gain energies several times 
of the plasma ram energy. This clearly shows that the 
accelerated particles are originated from the shock layer, 
not in the downstream region of the shock. This is con- 
sistent with previous one-d imensional hybrid simulations 
(|Kucharek fc Scholeil[T99l[ ). 

Since field lines themselves meander in space, it is dif- 
ficult to demonstrate that particles move off of field line 
simply by examining individual trajectories. To address 
this, in the upper panel of Figure [7] we show a scat- 
ter plot for particles with energies larger than 1.5Ei 
from the one-dimensional simulation at fldt = 100.0 
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Fig. 8. — The probability distribution function (PDFs) for the 
transverse displacements of the particles in the downstream region. 
The black solid line is for the test-particles originally along y = 
20c/uipi and z = 20c/ujpi in the three-dimensional simulation and 
the red dashed line is for the one-dimensional simulation. The 
PDFs have been normalized using the counted number of particles. 



in the x-z plane that were initiahy all on the same 
magnetic field line. Thi s m ethod has been used by 
I Giacalone k Jokipiil ()1994[) and lGiacaloni (|1994f ) to show 
that in one-dimensional and two-dimensional simula- 
tions, the motions of charged particles are restricted. It 
can be seen that in the one-dimensional simulation, the 
distribution of the particles is along a fluctuating field 
line, meaning that they arc restricted on the field line. In 
the three-dimensional simulation, we inject 120000 test- 
particles (the same number of particles as in the one- 
dimensional simulation) at the beginning of the simula- 
tion and track their motions during the simulation. The 
test-particles arc initially evenly distributed along line at 
y = 20c/ojpi and z — 20c/ujpi in the x direction (along 
the same field line). In the bottom panel of Figure [71 
we show the displacements of those particles with en- 
ergies larger than 1.5i?i at fldt = 120.0 similar to the 
one-dimensional case. We find that the distribution of 
the particles has no obvious struct ure, indicating they 
move off t heir origina l field line. iGiacalone fc Jokipiil 
(|1994f) and iGiacalonj ()1994f ) have shown that in one- 
dimensional and two-dimensional numerical simulations, 
the motions of charged particles are restricted on their 
original field lines. Here we show that this restriction is 
indeed removed in the three-dimensional hybrid simula- 
tion. We also observe that some test-particles can move 
a large distance (> 40c/ujpi) transverse to the magnetic 
field in the three-dimensional simulation. 

In Figure [8] we plot the probability distribution func- 
tions (PDFs) for the transverse displacement from the 

initial location ^^Ay^ + Az^ from the two cases shown in 
Figure [71 The black solid line represents the results from 
the three-dimensional simulation and the red dashed 
line represents the results from the one-dimensional sim- 
ulation. The PDFs have been normalized using the 
counted number of particles. It can be seen that a sig- 
nificant fraction of the particles can move a large dis- 
tance (> 4:0c/u!pi) in the transverse direction because 



Fig. 9. — Downstream energy spectra of protons in the three- 
dimensional simulation (black solid line) and the one-dimensional 
simulation (red dashed line) at Qcit = 120.0. 

the particles can move off their original field line. In the 
one-dimensional simulation, the transverse motion is re- 
stricted within about 20c/LUpi. In our case the field line 
meandering due to the large-amplitude magnetic fiuctua- 
tion can cause some particles to move transverse to the x 
direction. However, the results for the one-dimensional 
simulation (Figu re [71 and Figure [8]) are consistent with 
(|Giacalonelll994[) , who demonstrated the particles have a 
restricted cross-field motion in one- and two-dimensional 
fields by simulating a quasi-perpendicular shock with 
Obu = 45°. 

In Figure [HI we show a comparison of the down- 
stream energy spectra at ildt = 120.0 between the one- 
and three-dimensional simulations. The black solid line 
shows the spectrum for the three-dimensional simulation 
and the red dashed line represents the spectrum for the 
one-dimensional simulation. In both of the two cases, 
a fraction of thermal protons are accelerated to more 
than several times the plasma ram energy, Ep. Since in 
the one-dimensional simulation we have a limited num- 
ber of particles, the spectrum only extends to about 200 
■miV^Q/2. For the three-dimensional simulation, the par- 
ticles can be energized up to 600 miV^Q/2. In both of 
the cases, the density of the accelerated particles with 
energy larger than 20 miV^Q/2 is about 5% of the down- 
stream plasma density, indicating an efficient injection of 
energetic particles. 

4. DISCUSSION AND CONCLUSIONS 

We performed three-dimensional hybrid (kinetic ions 
and fiuid electrons) simulations to investigate the initial 
acceleration of thermal protons at parallel shocks. We 
examined the trajectories of many individual particles 
and the spatial distribution of protons that comprise the 
high energy tail resulting from the acceleration of parti- 
cles at the shock. We find that the charged particles gain 
energy at the shock layer similar to the one-dimensional 
simulation. The particle can ride on the shock front 
and gain a large amount of energy. The results con- 
firm previous hybrid simulations that the initial accel- 
eration of charged particles is right at the shock front 



Particle Acceleration at Parallel Shocks 



7 



(IQuestI [19881 IScholer fc Terasawal [T990t IScholed [l990l: 
iKucharek fc Scholeii 119911: iGiacalone et alJIlOOa . even 
in a three-dimensional electromagnetic field. The result 
clarifies the injection process for diffusive shock acceler- 
ation. In addition, we find that the guiding center of 
accelerated particle can move off their original field lines 
in the three-dimensional simulation. 

The "injection problem" is a long standing problem for 
diffusive shock acceleration. The charged particles are 
required to be energetic enough and efficiently interact 
with magnetic fluctuations. The results of our study bear 
on this problem. In this study we focus on the mecha- 
nism involved in accelerating thermal or low-energy par- 
ticles at the shock front. Previous studies have pro- 
posed two different mechanisms for the injection of low- 
energy ions at parallel s hocks, i.e., pa rticle energization 
by downstream heating (Ellison"198 lD and by reflection 
at the shock layer ([Quest 198&.Scholer fc Terasawal 19901: 
Scholei<[T990t IKucharek fc Scholedll991l: IGiacalone eTall 



199^. Although previous hybrid simulations have found 



that the initial energization is due to the ion reflection 
and acceleration at the shock layer, the results were ob- 



tained only for one-dimensional simulations and occa- 
sio nally two-dimensional simulations. As pointe d out 
bvlJokipii et al.l (I1993D . IGiacalone fc Jokipiil ()1994D . and 
iJones et al.l (|1998f ). in a magnetic field that has at least 
one ignorable coordinate, the motions of charged parti- 
cles are restricted on their original field lines of force. In 
this sense, the previous simulations are not conclusive. 
In this work, we have demonstrated that the reflection 
at the shock layer is the mechanism for accelerated parti- 
cles gaining their initial energy in a self-consistent, three- 
dimensional field. However, our results are not consistent 
with the "thermal leakage" model, which assumes that 
the thermal particles gain the first amount of energy by 
thcrmalization process in the downstream region. 
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